Collective excitation of Bose-Einstein condensates in the transition region between 3D 

and ID 
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We measure the frequency of the low m=0 quadrupolar excitation mode of weakly interacting 
Bose-Einstein condensates in the transition region from the 3D to the ID mean-field regime. Various 
effects shifting the frequency of the mode are discussed. In particular we take the dynamic coupling 
of the condensate with the thermal component at finite ternperature into account using a time- 
dependent Hartree-Fock-Bogoliubov treatment developed in U|. We show that the frequency rises 
in the transition from 3D to ID, in good agreement with the theoretical prediction Q. 
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One-dimensional (ID) quantum degenerate Bose gases 
have recently attracted considerable theoretical and ex- 
perimental interest Q. On the one hand this interest is 
due to their remarkable physical properties which are ab- 
sent in three-dimensional systems. On the other hand a 
rapid advance in trapping techniques for ultracold gases 
has put these systems within experimental reach. In par- 
ticular optical lattices Q, optical dipole traps and 
atom chips 6] have recently been used to realize such 
low-dimcnsional systems. For further experiments under 
these conditions a good understanding of the transition 
between the 3D and the ID regime is therefore of crucial 
importance. In this paper this transition region is char- 
acterized experimentally by monitoring the oscillation of 
a quantum degenerate Bose gas. 

From a fundamental point of view, one of the most 
striking features of ID quantum degenerate Bose gases 
is the predominant role played by quantum fluctuations. 
For spatially homogeneous ID systems, fluctuations of 
the phase rule out the existence of any off-diagonal long- 
range order (ODLRO) Q even at temperature T — 
The flnite size of trapped ID gases however al- 
lows for a rich variety of possible scenarios, including 
true phase-coherent Bose-Einstein condensates (BECs) 
as well as phase-fluctuating BECs, the so-called quasi- 
condensates 0. 

The behavior of these ID Bose gases is governed by 
the ratio of interaction and kinetic energy 
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where is the ID atomic density, m the atomic mass 
and giD the ID coupling constant. Since this ratio scales 
as 1/n these gases counter-intuitively become more non- 
ideal when the density is decreased. The fascinating 
features of these systems have led to continued theoret- 
ical interest over the past decades. For homogeneous 
ID Bose gases with short-range interactions, the ground 
state excitation spectrum and thermodynamic 
properties ^3 of the system can be determined with 



a Bethe Ansatz for arbitrary values of 7. In the case 
of trapped systems the equation of state can be found 
by combinin g th is approach with the local density ap- 
proximation [131 . For high densities the system is in 
the weakly interacting regime, where it can be well de- 
scribed in the frame of mean-field theories and ODLRO 
is present. On the contrary, for low densities the sys- 
tem enters the strongly interacting or strongly correlated 
regime, where a description by mean-field theories fails 
and ODLRO is strongly reduced. In the limit 7 00 the 
system is equivalent to a ID gas of impenetrable bosons, 
the so-called Tonks- Girardeau (TG) gas J^, where the 
bosonic particles effectively acquire fermionic properties. 

To experimentally prepare a ID gas in a harmonic trap 
with cylindrical symmetry one has to fulfill the condition 



(2) 



where u± denotes the radial trap frequency and ^ is the 
chemical potential of the ensemble. If this condition is 
fulfilled, neither the thermal nor the interaction energy 
is sufficient to affect the radial shape of the ground-state 
wave function. The radial degree of freedom is then 
frozen, since the atoms are confined to the ground state 
of the radial trapping potential. 

The first experimental realization of a weakly inter- 
acting ID BEG was achieved in a magnetic trap with 
very high aspect ratio and detected by a change in the 
ballistic expansion fl^l H^ . The tight radial confine- 
ment required to fulfill Q can also be provided by mag- 
netic microtraps IT?] or optical wave guides Such 
waveguides have allowed for the first realization of stable 
matter- wave bright solitons ^IQ]. Two-dimensional op- 
tical lattices 20] offer the possibility to realize arrays of 
ID gases which are coupled via tunneling between the lat- 
tice sites. In such a lattice gas a theoretically predicted 
reduction of the three-body recombination rate within 
the strongly correlated regime |2ll | has been experimen- 
tally observed [l^. Recently, seminal experiments have 
reached the TG regime in 2D optical lattices [2^ E3 | . 
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In this paper, we experimentally study collective ex- 
citations of weakly interacting BECs in the transition 
region between 3D and ID, i.e. for values of the chem- 
ical potential /i ~ hio±. The experimental investigation 
of collective excitations allows for quantitative tests of 
the underlying theoretical description. We focus on the 
observation of the low m = quadrupolar mode. Collec- 
tive modes were investi gate d experimentallv and theoret- 
ically in 3D BECs El lH El 113, EH and in ar- 
rays of ID gases confined in a 2D optical lattice |H E^ ■ 
Our experimental investigation in the transition region 
reveals clear signatures of one-dimensional behavior and 
is complementary to experiments on the phase-coherence 
properties of the system. 

The dynamics of the low m=0 mode consists of an 
out-of-phase oscillation of the axial and radial diameters 
of the cloud. For Bose gases confined in cylindrically 
symmetric traps the hydrodynamic equations of super- 
fluids 35] allow for the analytic determination of the 
mode frequency at T = 0: For the case of a very elon- 
gated but 3D weakly interacting BEC the frequency is 
u) = ■y/5/2 Wz js^, where ujz is the axial trap frequency. 
In the transition from the 3D mean-field regime to the 
ID mean-field regime, the radial dynamics of the con- 
densate freezes out, and accordingly the radial compo- 
nent of the low TO=0 mode vanishes. In this process, the 
mode frequency rises to uj — V^ujz jHls^. For the TG 
gas the frequency increases to uj — 2 ujz [Sg , identical to 
a non-interacting Fermi gas. The frequency in the in- 
termediate regimes has been numerically determined by 
using a combination of hydrodynamic equations with a 
sum-rule approach |3] and later by exclusively relying on 
hydrodynamic equations js^- Note that the analytic so- 
lutions given above depend on the axial trap frequency 
only, whereas the mode frequency in the intermediate 
regimes depends on the radial confinement as well 0. 

Our experiments were performed with ^^Rb Bose- 
Einstein condensates in the \F = l,mF — —1) hyperfine 
ground state confined in a strongly elongated magnetic 
trap. The axial trap frequency was 3.40 Hz and the radial 
trap frequency was varied between 265 Hz and 385 Hz, re- 
sulting in aspect ratios as low as A = uJz/uj± ~ 1/113. 
Further details of our experimental apparatus can be 
found elsewhere |4oj . In a trap with fixed geometry, the 
ID criterion reduces to a condition on the number 
of atoms in the sample. Neglecting the contribution of 
thermal atoms to the chemical potential, the gas is ID if 
the number of condensed atoms A^o is smaller than [l5| 



32 



225 a a i 



(3) 



where a± and denote the harmonic oscillator length in 
the radial and axial directions and a is the s-wave scat- 
tering length. For our trap parameters iVg^ « 4000. To 
approach this regime we reduce the number of atoms by 
adjusting the end of the forced evaporation ramp. We are 
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FIG. 1: Absorption images after 30 ms time of flight for vari- 
ous hold times r in the magnetic trap (indicated on the right) . 
The oscillation of the quadrupolar mode is clearly visible. The 
axial displacement of the cloud is due to the simultaneous ex- 
citation of the dipole mode. Axial density modulations ob- 
servable in the absorption images are related to the presence 
of phase fluctuations in the condensate [s^. We have verti- 
cally compressed the images for graphical representation. 



restricted by our detection system to a minimum number 
of A 2000 atoms. For this minimum number of atoms 
fi/huj± « 0.75, thus we obtain ensembles on the border 
of one-dimensionality. 

We performed the following experimental procedure: 
Laser-cooled atoms were loaded into a moderately elon- 
gated magnetic trap with A ~ 1/25. Radio- frequency 
evaporative cooling was performed to obtain tempera- 
tures just above the transition temperature. Then the 
trap was axially decompressed to reach the desired as- 
pect ratio, and the final evaporation ramp to obtain BEC 
was performed in the strongly elongated trap. The low 
m=0 mode was excited by modulating the current of the 
magnetic trap for five periods at a frequency close to the 
expected value of the mode frequency uj. The conden- 
sate was allowed to oscillate in the magnetic trap for a 
hold time r. Subsequently, the trapping potential was 
switched off and the atomic density distribution was de- 
tected by resonant absorption imaging after 30 ms time 
of flight. A bimodal fit to the density distribution was 
used to extract the aspect ratio, the total and condensed 
particle numbers and the temperature of the ensemble. 
By varying the hold time r in the magnetic trap the oscil- 
lation was stroboscopically monitored as shown in Fig.^ 
Since the mode frequency depends on the dimensional- 
ity of the system, the data was sorted according to the 
parameter P = NQXa/a±, which appropriately describes 
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FIG. 2: Aspect ratio of the clouds after time of flight as a 
function of the hold time in the magnetic trap. The fit to the 
data gives a mode frequency of w = 5.41 ± 0.02 Hz. The data 
points originate from clouds with an average dimensionality 
parameter P = 13.3 and reduced temperature T/Tc = 0.34. 

the degree of dimensionality for our purposes G|. Fig- 
ure 12 shows the aspect ratio of clouds within an interval 
of the parameter P versus the hold time r. To extract 
the frequency we fit this data with a damped sinusoidal 
function. 

The frequency rise which accompanies the transition 
from the 3D to the ID mean-field regime constitutes a 
shift of only 10%. Due to the small size of this effect 
other influences on the oscillation frequency have to be 
considered carefully. In principle, frequency shifts can be 
caused by effects beyond the mean-field approximation, 
by large oscillation amplitudes and by the finite temper- 
ature of the system. 

The frequency shift due to corrections to the mean-field 
approximation can be analyzed in the frame of the hydro- 
dynamic theory of superfluids. Using the first quantum 
correction to the Bogoliubov equation of state, this shift 
of the lowest quadrupolar mode at T = is |4l| 
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where /(A) is a function that depends on the geometry of 
the trap (see Eqn. (17) in 01). Since the gas parameter 
a^n(0) is very small for our experimental conditions, this 
frequency shift is below the 0.1% level and thus negligibly 
small. 

The hydrodynamic equations are also well suited to 
investigate the effect of large oscillation amplitudes at 
T = 0. Large amplitudes cause nonlinear coupling be- 
tween the normal modes of the system and shifts of the 
oscillation frequency. The frequency shift for the low 
171 = quadrupolar mode is 



where is the relative oscillation amplitude of the con- 
densate length in the trap and (5(A) is a factor depending 
on the trap geometry (see Eqn. (22) in ji^l)- The ampli- 
tude Az can be extracted from the oscillation amplitude 
in the time of flight images Aj*-*^ by using scaling the- 
ory For our parameters we obtain A^ ~ A™^/V2 
and the oscillation amplitudes in the trap are A^ < 20% 
for our measurements. Hence the corresponding fre- 
quency shift is limited to 6uj{Az)/lj < 0.5% and is thus 
clearly smaller than the expected shift due to dimensional 
effects. 

Let us now turn to the most important correction to 
the mode frequency, which is caused by the finite tem- 
perature of the system. Early experiments on collec- 
tive excitations in 3D systems investigated effects due 
to the finite temperature of the sample, including fre- 
quency shifts with respect to the zero temperature pre- 
dictions 0, 01 . On the theory side various approaches 
have been proposed to include finite temperature into the 
theory of weakly interacting trapped gases. To describe 
the observed features, such as the damping of excitations 
or the dynamic coupling between the condensed and ther- 
mal components, time-dependent mean-field schemes are 
appropriate. In this paper, we rely on the linearized time- 
dependent Hartree-Fock-Bogoliubov approach proposed 
in [l| . Within this approach the coupled equations for the 
dynamics of the condensate and the thermal component 
are not solved self-consistently, but perturbatively up to 
second order in the coupling constant g. Thus the non- 
physical gap in the self-consistent static quasi-particle 
excitation spectrum is replaced by a well-behaved gap- 
less Bogoliubov type of spectrum. Moreover at T = 
this approach recovers the first order quantum correc- 
tion and therefore includes the frequency shift due 
to corrections to the mean-field approximation. 
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FIG. 3: Normalized frequency difference between our mea- 
surements and the r=0 prediction in 01 as a function of the 
reduced temperature. The lines indicate the theoretical fre- 
quency corrections for an interaction parameter of rj = 0.33 
(solid) and rj = 0.41 (dotted) according to 
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P = No:^a/a, 

FIG. 4: Measurement of the mode frequency uj for the tran- 
sition between the 3D and the ID mean-field regimes. The 
frequency shifts due to finite temperature and large oscillation 
amplitudes are included. The solid line shows the theoretical 
prediction 21- The dashed vertical line indicates the transi- 
tion point between the two regimes according to Q. 



Figure 13 shows a comparison between the measured 
frequency shifts with respect to the T = calculation 
in [3 and the theoretical prediction for the shifts due 
to finite temperature according to Q. Tlie calculated 
frequency shift depends on temperature and on the pa- 
rameter J] = [ij (kBTc) which describes the interactions 
in a 3D Bose gas 1]. Here Tc denotes the critical tem- 
perature for a non-interacting gas. Due to the variation 
of the atom number in the experiments, the correspond- 
ing values for rj range from 0.33 to 0.41. The theoretical 
prediction for these cases is shown in Fig. |2| Despite 
the low atom numbers used in some of the measurements 
this theory, valid in the thermodynamic limit, shows good 
agreement with the measured frequency shifts. Since the 
finite-temperature shifts are of the order of a few percent, 
they constitute an important contribution to the mode 
frequency. For a subsequent correction of the mode fre- 
quency, the parameter rj and the resulting frequency shift 
was calculated for each data point. 

To compare our finite-temperature measurements with 
the zero-temperature prediction in we correct our 
data for the finite-temperature shift and the smaller shift 
due to the oscillation amplitude. Figure0]shows the mea- 
sured mode frequencies, including these corrections. The 
variation in the number of particles allowed for measure- 
ments covering almost the entire transition region be- 
tween the 3D and ID mean-field regimes. Within the 
experimental error, our data shows good agreement with 
the predicted frequency dependence at zero temperature. 
The error increases towards small values of the dimen- 
sionality parameter P due to the small particle numbers 
involved in these measurements. A small discrepancy 
between the theoretical prediction and the measurement 



might be due to the loss of atoms during the oscillation 
or systematic uncertainties in the determination of the 
number of atoms. 

In conclusion, we have measured the increase of the 
oscillation frequency of the low m—0 quadrupolar mode 
in the crossover region between the 3D and ID regimes. 
To compare our measurements with the zero-temperature 
prediction, various effects that can lead to shifts of this 
frequency were carefully evaluated. We have identified 
finite-temperature effects as an important contribution 
to the frequency shift. Including frequency corrections 
due to the finite temperature and the oscillation ampli- 
tude, our data shows good agreement with the theoretical 
prediction. The observed frequency increase constitutes a 
clear signature of the onset of one-dimensionality. These 
results confirm that mode frequency measurements pro- 
vide a sensitive probe of the dimensionality of quantum 
degenerate Bose gases. 
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cussion with S. Giorgini and M. Lewenstein. This work 
was supported by the Deutsche Forschungsgemeinschaft 
within the SFB 407. 



* schulteQiqo.uni-hannover.de 
[1] S. Giorgini, Phys. Rev. A 61, 063615 (2000). 
[2] C. Menotti and S. Stringari, Phys. Rev. A 66, 043610 
(2002). 

[3] L. Pricoupenko, H. Perrin and M. Olshanii, Eds., J. Phys. 

IV, 116 (2004) 
[4] I. Bloch, J. Phys. B 38, 629 (2005). 

[5] R. Grimm, M. Weidemiiller, and Yu.B. Ovchinnikov, 

Adv. At. Mol. Opt. Phys. 42, 95 (2000). 
[6] J. Fortagh, S. Kraft, A. Giinther, C. Triick, P. Wicke, 

and C. Zimmermann, Opt. Comm. 243, 45 (2004). 
[7] J. W. Kane and L. P. Kadanoff, Phys. Rev. 155, 80 

(1967); J. Math. Phys. (N.Y.) 6, 1902 (1965); L. Reatto 

and G. V. Chester, Phys. Rev. 155, 88 (1967). 
[8] M. Schwartz, Phys. Rev. B 15, 1399 (1977); F. D. M. 

Haldane, Phys. Rev. Lett. 47, 1840 (1981). 
[9] D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, 

Phys. Rev. Lett. 85, 3745 (2000). 
[10] E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 
[11] E. H. Lieb, Phys. Rev. 130, 1616 (1963). 
[12] C. N. Yang and C. P. Yang, J. Math. Phys. (N. Y.) 10, 

1115 (1969). 

[13] V. Dunjko, V. Lorent, and M. Olshanii, Phys. Rev. Lett. 

86, 5413 (2001). 

[14] M. Girardeau, J. Math. Phys. (N. Y.) 1, 516 (1960). 
[15] A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, 

T. L. Gustavson, J. R. Abo-Shaeer, A. P. Chikkatur, S. 

Gupta, S. Inouye, T. Rosenband, and W. Ketterle, Phys. 

Rev. Lett. 87, 130402 (2001). 
[16] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, T. 

Bourdel, J. CubizoUes, and C. Salomon, Phys. Rev. Lett. 

87, 080403 (2001). 

[17] W. Hansel, P. Hommelhoff, T. W. Hiinsch, and J. Reichel, 



5 



Nature 413, 498 (2001); J. Fortagh, H. Ott, S. Kraft, A. 
Giinthcr, C. Zimmcrmann, Appl. Phys. B 76, 157 (2003). 
[18] K. Bongs, S. Burger, S. Dettrner, D. HcUweg, J. Arlt, W. 
Ertmer, and K. Sengstock, Phys. Rev. A 63, 031602 (R) 
(2001). 

[19] L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. Cu- 

bizolles, L. D. Carr, Y. Castin, and C. Salomon, Science 

296, 1290 (2002); K. E. Strccker, G. B. Partridge, A. 

G. Truscott, and R. G. Hulet, Nature 417, 150 (2002); 

B. Eierrnann, Th. Anker, M. Albiez, M. Taglieber, P. 

Treutlein, K.-P. Marzlin, and M. K. Oberthaler, Phys. 

Rev. Lett. 92, 230401 (2004). 
[20] M. Greiner, I. Bloch, O. Mandel, T. W. Hansch, and T. 

Esshnger, Phys. Rev. Lett. 87, 160405-1 (2001). 
[21] D. M. Gangardt and G. V. Shlyapnikov, Phys. Rev. Lett. 

90, 010401-1 (2003). 
[22] B. L. Tolra, K. M. O'Hara, J. H. Huckans, W. D. PhiUips, 

S. L. Rolston, and J. V. Porto, Phys. Rev. Lett. 92, 

190401 (2004). 

[23] B. Paredes, A. Widera, V. Murg, O. Mandel, S. Filing, 
L Cirac, G. V. Shlyapnikov, T. W. Hansch, and I. Bloch, 
Nature 429, 277 (2004). 

[24] T. Kinoshita, T. Wenger, D.S. Weiss, Science 305, 1125 
(2004). 

[25] D. S. Jin, J. R. Ensher, M. R. Matthews, C. E. Wieman, 

and E. A. Cornell, Phys. Rev. Lett. 77, 420 (1996). 
[26] M.-O. Mewcs, M. R. Andrews, N. J. van Druten, D. M. 

Kurn, D. S. Durfee, C. G. Townscnd, and W. Ketterle, 

Phys. Rev. Lett. 77, 988 (1996). 
[27] C. Fort, M. Prevedelh, F. Minardi, F. S. Catahotti, L. 

Ricci, G. M. Tino, and M. Inguscio, Europhys. Lett. 49, 

8 (2000). 

[28] F. Chevy, V. Bretin, P. Rosenbusch, K. W. Madison, and 
J. Dalibard, Phys. Rev. Lett. 88, 250402 (2002). 

[29] R. Onofrio, D. S. Durfee, C. Raman, M. Kohl, C. E. 
Kuklewicz, and W. Ketterle, Phys. Rev. Lett. 84, 810 
(2000). 



[30] S. Stringari, Phys. Rev. Lett. 77, 2360 (1996). 

[31] M. Edwards, P.A. Ruprecht, K. Burnett, R.J. Dodd, and 

Charles W. Clark, Phys. Rev. Lett. 77, 1671 (1996). 
[32] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 

Rev. Lett. 91, 250402 (2003). 
[33] P. Pedri, and L. Santos, Phys. Rev. Lett. 91, 110401-1 

(2003). 

[34] S. Dcttmcr, D. HcUweg, P. Ryytty., J. J. Arlt, W. Ert- 
mer, K. Sengstock, D. S. Pctrov, G. V. Shlyapnikov, H. 
Kreutzmann, L. Santos, M. Lcwenstcin, Phys. Rev. Lett. 
87, 160406 (2001); S. Richard, F. Gerbier, J. H. Thywis- 
sen, M. Hugbart, P. Bouyer, A. Aspect, Phys. Rev. Lett. 
91, 010405 (2003). 

[35] P. Nozieres and D. Pines, The Theory of Quantum Liq- 
uids, (Addison- Wesley, Reading, MA, 1990), Vol. H. 

[36] S. Stringari, Phys. Rev. A 58, 2385 (1998). 

[37] Tin-Lun Ho and M. Ma, J. Low Temp. Phys. 115, 61 
(1999). 

[38] A. Minguzzi, P. Vignolo, M. L. Chiofalo, and M. P. Tosi, 

Phys. Rev. A 64, 033605 (2001). 
[39] J. N. Fuchs, X. Lcyronas and R. Combescot, Phys. Rev. 

A 68, 043610-1 (2003). 
[40] D. Hellweg, L. Cacciapuoti, M. Kottke, T. Schulte, K. 

Sengstock, W. Ertmer, and J. J. Arlt, Phys. Rev. Lett. 

91, 010406 (2003). 
[41] L. Pitaevskii and S. Stringari, Phys. Rev. Lett. 81, 4541 

(1998). 

[42] F. Dalfovo, C. Minniti, and L. P. Pitaevskii, Phys. Rev. 

A 56, 4855 (1997). 
[43] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996). 
[44] D. S. Jin, M. R. Matthews, J. R. Ensher, C. E. Wieman, 

and E. A. Cornell, Phys. Rev. Lett. 78, 764 (1997). 
[45] D. M. Stamper-Kurn, H.-J. Miesner, S. Inouye, M. R. 

Andrews, and W. Ketterle, Phys. Rev. Lett. 81, 500 

(1998). 



